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Abstract 

We perform a systematic study of the phenomenology associated to models where the 
dark matter consists in the neutral component of a scalar SU(2)l n-uplet, up to n = 7. 
If one includes only the pure gauge induced annihilation cross-sections it is known that 
such particles provide good dark matter candidates, leading to the observed dark matter 
relic abundance for a particular value of their mass around the TeV scale. We show 
that these values actually become ranges of values - which we determine - if one takes 
O |' into account the annihilations induced by the various scalar couplings appearing in these 

models. This leads to predictions for both direct and indirect detection signatures as 
a function of the dark matter mass within these ranges. Both can be largely enhanced 
by the quartic coupling contributions. We also explain how, if one adds right-handed 
neutrinos to the scalar doublet case, the results of this analysis allow to have altogether a 
viable dark matter candidate, successful generation of neutrino masses, and leptogencsis 
in a particularly minimal way with all new physics at the TeV scale. 



o 

^ : 1 Introduction 

O , 

Q\ , There are many possible dark matter (DM) candidates and a systematic study of all possi- 

bilities, and associated phenomenology, is not conceivable. However if one takes as criteria 
the minimality of the model, in terms of the number of new fields and parameters, such a 
systematic study becomes feasible. Such approach is different and complementary to the 
ones that led to theories with e.g. Supersymmetry [1-3] or Universal Extra Dimensions [4-6], 
which were invented as an attempt to address and solve other fundamental questions such as 
the Hierarchy problem, and where the number of parameters and possibilities can be huge. 
Another criterion of selection one can consider is the predictivity and the testability of the 
model in current and future accelerators, and direct or indirect DM detection experiments. 
Particularly simple possibilities along these lines of thought arise if one adds to the Standard 
Model only one extra SU(2)l singlet or multiplet, scalar or fermion, containing a neutral 
DM candidate field. The stability of the DM is usually achieved in this case by introducing 
a Z2 parity symmetry, under which the extra multiplet is odd and all the SM particles are 
even. 

Several possibilities of this kind, such as the scalar singlet [7-14], the fermion singlet 
[15,16], the scalar doublet (in the "Inert Doublet Model" [13,17-24]), the fermion doublet 
candidate [25,26], etc, have already been explored and they offer a rich phenomenology. A 
systematic study has been performed in Ref. [25] for any multiplet from the doublet up to the 
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7-plet. Multiplets offer the advantage that they could be potentially produced at colliders 
through gauge interactions. In this analysis the relic density of such DM candidate has been 
calculated considering all annihilation processes induced by the known SU{2)l x U(l)y gauge 
interactions. This framework is particularly predictive, the only free parameter is the DM 
mass, rriDMi an d the observed relic density can be obtained for only one value of this mass. 

Considering only the gauge induced processes in such a way is fully justified for an extra 
fermion multiplet because no other renormalizable interaction with the SM particles can be 
written. However, for a scalar multiplet this assumption is not at all automatic as quartic 
scalar interactions involving both the scalar multiplet and the Brout-Englert-Higgs doublet 
are perfectly allowed. Therefore, the analysis of Ref. [25] for scalar multiplets does not hold 
if these scalar couplings are not suppressed. 

In this paper, we study in a systematic way the rich phenomenology which arise if one 
includes the effects of the quartic couplings for all scalar multiplets up to the 7-plet. This 
we do in the high mass regime, that is to say for tudm > m w> where the observed relic 
density is obtained for annihilation cross-sections oc l/m 2 DM (which is typical of the large 
DM mass asymptotic regime). We show in particular that, due to a large enhancement of the 
(co) annihilation of DM into gauge bosons driven by the scalar couplings, the latter cannot 
be ignored unless they are much smaller than the gauge couplings. Moreover, due to these 
quartic couplings, and without fine-tuning, a large range of values of uidm is compatible with 
the observed DM relic abundance. These contributions also enhances the predicted fluxes for 
direct and indirect detection searches. 

The case where the multiplet is a doublet, known as the Inert Doublet Model (IDM) 
has already been extensively studied in the literature. A detailed analysis of the high mass 
regime was however missing and we provide it here. The phenomenology of this model is 
particularly rich because it depends on the interplay of three different scalar quartic couplings. 
The phenomenology of the higher multiplet case on the other hand in fine depends on only 
one quartic coupling, A3, which renders these cases particularly constrained and predictive. 
In this work, only higher multiplet models allowed by current direct detection constraints 
will be considered. This limits us to odd dimension n-uplets with zero hypercharge and 
n = 3, 5 and 7. For these models the high mass regime, moj/ > niw, which we study is the 
only possible one. 

We also present in this paper an intriguing possible consequence of our results for the 
doublet model: in agreement with the DM constraints, if, in order to explain the neutrino 
masses, one adds to this model right-handed neutrinos, it is possible to induce in a particularly 
simple way baryogenesis through leptogenesis with all new physics around the TeV scale. 

The paper is organized as follows. We first present the inert doublet and the higher mul- 
tiplet models in Section [2j Aside from fixing notations and definitions, a discussion on the 
number of relevant quartic couplings for the higher multiplets is made. Predictions for the 
relic density are made in Section [3l both numerically and (in the instantaneous freeze-out 
approximation) analytically. The latter method allows to show the enhancement of the scalar 
coupling contribution in the various cross-sections, in particular the important coannihilation 
ones. For the doublet case, maximal mass splittings between the DM doublet components 
compatible with the WMAP constraint are given as a function of the DM mass. For higher 
multiplets, this constraint fixes the value of A3 as a function of the DM mass. A discussion is 
also made about the consequences of having, for very heavy DM candidates, freeze-out before 
the electroweak phase transition. In Section U predictions on the DM-nucleon elastic scat- 
tering cross-section relevant for direct detection searches are made, and compared to current 
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experimental limits and projected reaches of future experiments. In Section predictions for 
various indirect detection signals are discussed. The possibility of resonances [27, 28] is reex- 
amined in light of the enlarged mass range of the DM candidate. Photon and neutrino fluxes 
from the galactic center are compared with the sensitivity of current telescopes (FERMI and 
KM3net). The fluxes of charged antimatter cosmic rays (positrons and antiprotons) are cal- 
culated with DarkSUSY and confronted with data. Finally, the extension of the doublet by 
right-handed neutrinos and consequences for neutrino masses and leptogenesis, are discussed 
in Section [6l Conclusions are drawn in Section [7[ Appendix A contains precise discussions, 
for the higher multiplet cases, on the most general scalar potential and the differences and 
the similarities between complex and real multiplets. Appendix B gives the complete set of 
(co) annihilation Feynman diagrams for all the models studied in this paper. 

2 Models 

2.1 Inert Doublet Model 

The Inert Doublet Model (IDM) is a two Higgs doublet model with a Z 2 symmetry. They are 
denoted by Hi and H%, Hi being the usual Brout-Englert-Higgs doublet. All SM particles 
are even under the Z 2 symmetry, while H2 is odd. This ensures the stability of the lightest 
member of H2, which will be the DM candidate, and prevents from flavor changing neutral 
currents (FCNC) [17]. We will assume that Z2 is not spontaneously broken, in particular, 
H2 does not develop a vacuum expectation value. In order to have a neutral component, 
the hypercharge of a scalar doublet is necessarily Y = ±1 (we choose to write the electric 
charge Q = T3 + Y/2). We conventionally assign +1 to the hypercharge of H2: one can 
write H2 = (H + (Hq + L4q)/\/2) t , similarly to the ordinary Higgs doublet, where Hi = 
(h+ ( V0 + h + iG )/V2) T . 

The most general renormalizable scalar potential with two doublets is given bjo 



V{Hi,H 2 ) = ^ 2 |#i| 2 + fJ- 2 2 \H 2 \ 2 + Ai|#i| 4 + A 2 |# 2 | 4 

+ A 3 |#i| 2 |# 2 | 2 + A 4 |if 1 t ^ 2 | 2 + y [{H\H 2 f + h.c. 



(1) 



After the electroweak symmetry breaking, Hi develops its vev, vq = — /xf/Ai — 246 GeV, and 
the scalar potential in the unitary gauge then becomes, 

V = -m 2 h h 2 + Xiv h 3 + ^Ai/1 4 

+i {\ Ho H 2 + X Ao A 2 + 2X Hc H + H~) (2v h + h 2 ) 

+ 1\2{H 2 + A 2 + 2H + H-) 2 , (2) 



1 For the doublet case, the introduction of the term {H\TiH\){H\TiH2) , where the r< are the SU(2) gener- 
ators, is redundant since 



(H\Hi){HlH2) + {HlnH^HlnHi) = 2\HlR. 



21 2 
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with a mass spectrum given by, 



m l = 2X i v o > 

m H = A + X H Vq , 

m A = lA + ^A v 2 o , 

m H+ = lA + Ah c ^o • (3) 

We have defined \h c = A3/2 and Xh ,a = (A3 + A4 ± As)/2. We will consider ifo to be the 
DM candidate (i.e. A5 < 0) though the results would be exactly the same for Aq changing 
the sign of A5. 

Some theoretical constraints first apply on these quartic couplings. To ensure that the 
potential is bounded from below, the vacuum stability (at tree-level) requires that, 

Ai, 2 > , 



Afro > A^o , \h c > - v A1A2 . (4) 

Other constraints are also imposed from past accelerator measurements. Indeed the extended 
scalar sector could bring corrections to electroweak precision test observables (EWPT). In 
particular, the variable T, which is a measure of the radiative corrections to myy /{jnz cos 0w), 
has been calculated for the IDM [19] 

AT ~ V2ix 2 av 2 ( mjJ+ ~ mAo ^ mH+ ~ mH ^ ■ ( 5 ) 

However, in the high mass regime considered in this work, mass splittings turn out to be 
small (see Section [3]), so that this constraint is not limiting. We find AT < 6.0 • 10 -3 , well 
below current experimental bounds [29,30]. 

The IDM has already been extensively studied in the literature. It has been shown that a 
viable DM candidate with the correct relic abundance can be obtained in three regimes, low- 
mass (mH Q <§C fnw) [13,23], middle-mass (mH — n^w) [19,21] and high-mass (rriH ^> m\y) 
[21,25]. Direct and indirect detection constraints were investigated in Refs. [13,19-22,31-33] 
and confrontation of the IDM in the muo < m w regime to colliders data and future prospects 
was done in [34,35]. In this paper, we provide a more detailed analysis of the high mass regime 
and show that the scalar coupling contribution can easily dominate over the gauge one and 
without fine-tuning (as results in Ref. [21] suggested). 

2.2 Higher Multiplet Models 

The procedure followed for the doublet above can be generalized for a multiplet of higher 
dimension. Let H n denotes this scalar multiplet, with n being the dimension of its represen- 
tation under SU{2)i. The relevant lagrangian for any of these objects coupled to the usual 
Higgs doublet H\ can be written as 

C = (D u H n y {D»H n )-V(H n ,H{) , (6) 

with the covariant derivative given by 

D^d.-igr^W^-igy^B^ . (7) 
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stands for SU(2)l generators in the representation n and Y is the hypercharge of H n . 
For H n to contain a neutral component H n °\ the hypercharge Y has to be odd (even) when 
n is even (odd). The most general renormalizable potential for H n is given by 



V(H n , Hx) = Vx{H x ) + fi 2 HlH n + ^ (h^H^ + A 3 (h\Hx) (h^H, 

A 



(8) 



where a sum over a is implicit in the last two terms. In this potential, A2,3 are equivalent 
to A2,3 in the doublet case while the A5 operator is equivalent to the sum of the A4 and A3 
operators in the doublet case. The A4 operator reduces to the A2 operator in the doublet 
case. It is important to notice that the A5 operator in the doublet case, which is responsible 
for the mass splitting between Hq and Aq, has no equivalent in higher multiplet dimension. 

Unlike the doublet case, the potential of Eq. ([8]) cannot give rise to a mass splitting 
between the real and the imaginary part of the neutral component of the multiplet, when 
this field is complex. IF Y ^ this would lead to a DM candidate with unsuppressed vector 
interactions with the Z boson, which is ruled out by direct detection limits. Unless some 
mechanism is advocated to create this mass splitting, this restricts the viable models to odd 
dimension multiplets with Y = (the coupling to Z is proportional to (T3 — Qsin 2 6w), with 
Q = T3 + Y/2). Notice that there are still two cases with Y = 0, depending on whether the 
multiplet is real or complex. The perturbativity of SU(2)l up to the Planck scale imposes 
n < 8 [25]. Therefore, the only possibilities are n = 3, 5, 7. 

As in the doublet case, for n = 3 and n = 5, a Z2 symmetry is necessary to ensure the 
stability of the DM candidate. In the case n = 7, this parity is unnecessary because the 
candidate is automatically stable. Indeed, no renormalizable or dimension 5 operator can be 
constructed to induce its decay into SM particles [25] . Moreover, an operator of dimension 6 
or higher would induce a lifetime of the order of the age of the universe or larger if the cutoff 
scale is set to the GUT scale. The DM candidate for n = 7 is accidentally stable, like the 
proton in the SM. 

Let us first analyze the case of the real multiplet models. In a suitable basis, as detailed 
in Appendix lA.il a real multiplet H n is written as 



/ \ 

a'(°) 



(9) 



where j n = (n — l)/2, Q in A^ corresponds to the electric charge, and A^ ^ = (A^))*. 
For real multiplets, the expression ( H n T^ H n ) is identically zero. Therefore the terms with 



coefficient A4 and A5 disappear from the potential Eq. (jSJ). As a consequence, there is only 
one scalar quartic coupling (A3) connecting H n to Hx- 

After the electroweak phase transition, the SM Higgs field develops its vev, (Hx) = vo/V^, 
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and the scalar potential in the unitary gauge becomes 
V = \rn\h 2 + Ait-o/i 3 + ^Ai/i 4 

+ I m 2 A (0)2 + £ m 2 A (Q) A (-Q) 

0<Q<jn 



\ 0<Q<j„ / 

\ 0<Q<j„ / 



(2t> /i + /i 2 ) 



+ i - AW a + \ • A iw; A i-«j 1 . (10 ) 



At tree-level, all the multiplet components have the same mass 

m§ = „ a + ^ . (11) 

At one-loop however, a mass splitting is generated by the coupling to gauge bosons and the 
charged components become slightly heavier than the neutral one [25], 



A«))-m(Al°')=Q 2 AM 9 , (12) 



m 
where 

AM g = gM w sin 2 ^ ~ (166 ± 1) MeV . (13) 

Notice that the scalar couplings of the potential do not modify these splittings, because they 
are identical for all the charged and the neutral components. As these one-loop splittings 
are small, it is a good approximation to consider all DM states as degenerate. Finally, the 
vacuum stability is ensured by the condition 

Ai, 2 > , 

A 3 > -V2AiA 2 . (14) 



The case of complex multiplets is analyzed in details in Appendix IA.21 It appears that 
the associated phenomenology is close to the real case. Without the introduction of some 
new symmetry U (1) under which H n is charged, a complex multiplet can be decomposed into 
two interacting real multiplets. 

In the presence of such a symmetry, all the degrees of freedom are degenerate at tree level 
except for the A5 term of Eq. (JSj) which induces an extra mass splitting (see Eq. ([71]) ) and 
lowers the mass of half of the charged components of H n . The neutral DM field stays the 
lightest only if A 5 < 2.2 • 1(T 2 . In the latter case, the model is similar to a real multiplet 
model, except for the doubling of the number of fields. This, in turn, reduces the threshold 
mass imposed by the relic density constraint by a factor y/2 which implies that scalar DM 
candidates lighter than the real case analyzed in what follows are still allowed. 



3 Relic abundance in the high-mass regime 



In this section, we show that the SU(2)l scalar multiplet extension of the SM can natu- 
rally lead to a multi-TeV DM candidate with the correct relic density. In this high-mass 
regime, coannihilations play a significant role. We will therefore start by briefly reviewing 
the formalism used to calculate the relic density. 



3.1 Freeze-out equations 

To calculate the DM relic abundance, we solve the Boltzmann equation for the total density 
of all the coannihilating species n = rii (we take i = for the lightest DM candidate and 
i > for the other species), 

^ + 3Hn = _ {(7effv){n 2_ n e q 2 ) ^ (15) 



where the effective thermal cross-section, given by 

M = £<^>^ (i6) 

is an average of the various thermal (co) annihilation cross-sections (a^v), weighted by equi- 
librium densities 

n^ = {m l T/2Trfl 2 e- m ^ T . (17) 

Although a full integration of this Boltzmann equation is needed in the general case to 
compute the relic density, an instructive and reliable estimate is derived from the so-called 
instantaneous freeze-out approximation, when poles or thresholds don't appear in the cross- 
sections [36]. For cold DM, we can develop the cross-sections in the non-relativistic limit, if 
av = A + Bv 2 , the corresponding thermal average is given by [37,38] 

(av) = a + b(v 2 ) = A + 6 ( B ~ : j]^ » ( 18 ) 

with x = mo/T. Then, the relic density is simply obtained as 

n , o 1-07 10 9 GeV" 1 . . 

VvMh 2 ~ m . (19) 

J{xF)gJ mp\ 
The post freeze-out annihilation integral J is given by 

J{XF ) = r ^s^dx , (20) 

and the freeze-out point xf is found by solving the equation 

0.0038 mpi g eff m (a eff v) 
xp = in — , (21) 

{g*x F yi 2 

where g e g = X^j / n o l * s the effective number of degrees of freedom. Usually, xf — 25. 
When all the DM species are degenerate, the equilibrium densities for all the states are 
equal, and the effective thermal cross-section {cr e ffv) is simply the average of (a^v) over all 
(co) annihilation channels. 
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3.2 Inert Doublet Model 



In a first step, we will consider the case without quartic couplings between H\ and H2, and 
derive the relic abundance constraint on the DM mass. Then, in a second step, we will 
show how the conclusions drawn in the step one are dramatically changed when these scalar 
couplings are present. 



3.2.1 IDM in the pure gauge limit 

When all the quartic couplings between H\ and H2 vanish (except A5 which we take tiny 
but non vanishing to avoid the direct detection problem above), all states are degenerate at 
tree-level. At one-loop, the neutral states remain exactly degenerate due to the Peccei-Quinn 
symmetry, up to the very small A5 contribution, and a splitting m H + — rriH — 350 MeV is 
induced between the charged and the neutral states [25]. Because of the smallness of these 
splittings for the annihilation cross-section, it is a very good approximation to consider all 
states as exactly degenerate and all quartic couplings as vanishing. In this limit, the DM 
species (co)annihilate into either known gauge bosons or fermions through an intermediate 
gauge boson. The corresponding Feynman diagrams are shown in Fig. [10] and Fig. [12] of 
the Appendix [B] The only free parameter is the DM mass, so that there is a one-to-one 
correspondence between the DM mass and the relic density. 

Let <To — Aq/v + Bq v be the effective cross-section in this limit, it is given by the average 
of the cross-sections for all the annihilation and coannihilation processes relevant for the relic 
density calculation. We obtain, at leading order, 

(3-2 4),°j (22) 

where 02 = <? 2 /47r is the weak coupling constant, s w = sin8]y, and c w = cos0\y- The zero 
velocity term Aq agrees with the result for n = 2 of Eq. (12) in Ref. [25], up to a factor one 
The velocity dependent term is mainly due to coannihilations, its coefficient -Bo is of 
the same order of magnitude as Aq. The analytical expression for Bq will not be given here, 
but we took it into account in numerical evaluations. 

It is interesting to notice that all the (co) annihilation cross-sections fall required 
by unitarity constraints. For annihilations into gauge bosons, this behavior is achieved after 
the cancellation of various diagrams whose amplitudes are connected by gauge invariance. Let 
us for instance examine in more details the process HqHq — ► ZZ. Naively, the contribution 
from the longitudinal modes Zl to the amplitude is enhanced by a factor m 2 Ho /m 2 z compared 
to the contribution from the transverse modes Zt- This would lead to an unacceptable 
behavior of the cross-section, a ~ m 2 H . Actually, a cancellation of the longitudinal parts 
occurs between the t and the u-channels on one hand, and the point-like interaction diagram 
("p"— channel) on the other hand. Notice that the t and the u-channels involve the propagator 
of Aq. When quartic couplings vanish, all DM states are degenerate, so that this cancellation 
is almost exact in the sense that the residual amplitude to longitudinal modes is given by 

n 2 Tn 2 

Ml = M(H H ->Z l Z l )~ / I . (23) 

4:C w m Ho 



2 The extra factor 1/2 in Ref. [25] could be due to a convention in the definition of the thermal average. 
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Figure 1: Dark matter relic abundance in the pure gauge limit as a function of the DM 
mass. Dashed (Solid) curve : Instantaneous freeze- out approximation without (with) velocity- 
dependent terms in o~v. Points : Output from MicrOMEGAs 



For the two transverse modes, there is no cancellation, their amplitude amounts to 

a 2 

M Ti = M(H H - Z Ti Z Ti ) ~ f_ (i = l,2) . (24) 

In the high mass regime, Ml "C M-Ti- For e.g. rriH ~ 550 GeV, Mh/Mn — 1.4 %. 
As a(H H — > ZZ) oc Sj|A^Ti| 2 + I-^lI 2 ) we see that the residual longitudinal amplitude 
is totally negligible. Therefore, in the pure gauge limit, gauge bosons produced by the 
annihilations of Hq are almost purely transverse. 

In Fig. [H we have plotted the DM relic density as a function of mass, assuming zero 
quartic couplings between Hi and Hi- The solid and the dashed curves correspond to the 
instantaneous freeze-out approximation with and without velocity dependent terms in (cry). 
As can be seen, these terms shift down the value of Qfjm^ 2 by only ~ 4%. Also shown 
are more exact points from a full integration of the Boltzmann equation, obtained with 
the MicrOMEGAs program [39]. The latest 5-years WMAP results, combined with baryon 
acoustic oscillations and supernovae data yield Odm^ 2 = 0.1131 ± 0.0034 [40]. We see that, 
in the absence of quartic couplings, the DM mass is determined by the relic density, 

m Ho = 534 ± 8.5 GeV (la) , (25) 

in agreement with the results of Ref. [25], up to the update of Odm^- 2 - At 3a the DM mass 
cannot be lighter than 508 GeV. It is worth noticing that the value of m# is quite sensitive 
to the precision at which Odm^ 2 is determined. Also, the approximate solid curve from Fig.[T] 
gives a slightly higher mass range, mu = 553 ± 8.5 GeV. The discrepancy is attributable to 
the instantaneous freeze-out approximation rather than to the values of the cross-sections in 
Eq. (SI. 

3.2.2 Effect of the quartic couplings 

When the scalar quartic couplings between H\ and H2 are switched on, the cross-section is 
affected in two ways. First, non-zero mass splittings between members of the inert doublet, 
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Eq. ([3]) will modify the amplitude of pure gauge diagrams of both annihilation and coanni- 
hilation processes. Second, a series of new annihilation and coannihilation processes which 
involve the usual Higgs particle appear, see Fig. [IT] of the Appendix IB1 

It is instructive to analyze how cross-sections grow with the quartic couplings. As a 
generic example, let us consider again the process HqHq — > ZZ. When a non zero mass 
splitting between Hq and Aq exists, the amplitudes to longitudinal modes from the t, u and 
p-channels do not cancel exactly anymore and is not suppressed anymore by a m 2 z /m 2 Ho as 
in Eq. (|23p . Instead, a contribution proportional to {m\ — m^) remains (in the high-mass 
regime squared mass splittings are small with respect to m 2 H o but not necessarily with respect 
to m%, as we will see), 

1 1 1 

. n z mir m An — mi, 
24 m\ m 2 Ho 

Furthermore, there is also a new contribution from the Higgs exchange in the s-channel 
HqHq — > h* — ► ZZ. The amplitude of the s-channel is proportional to \ho v q = {mj lQ — fj, 2 ). 
When added to the t, u and p-channels, a further cancellation takes places, and the total 
amplitude to longitudinal modes due to quartic couplings is finally given by 

M X „ 9 2 m H {(^Ao-^o) + (^q-^)} 
L 2c l rn% m H 

= 2A Ao (27) 

to be compared with Eqs. (|23M24p . Corrections to transverse modes are negligible, because 
they are smaller by a factor m? z /m 2 HQ . Therefore, gauge bosons produced by the scalar 
quartic couplings are almost purely longitudinal in the high mass regime, whereas those from 
the gauge interactions alone are almost purely transverse. As a consequence, the annihilation 
cross-section can only grow (if we neglect the tiny residual amplitude of Eq. (123 j) ) when 
scalar quartic couplings are switched on. The scalar coupling contribution to the cross- 
section <j(HqHq — > ZZ) becomes comparable to the gauge one for \a — g 2 /(2\/2c 2 u ) ~ 0.2. 
This corresponds to a small value of the splitting \mA — M2I ~ m w /A*2- 

The above analysis for the process HqHq — > ZZ serves as a demonstration that the various 
(co) annihilations cross-sections can only grow with the splittings between fj,2, Hq, Aq and H c 
(H c stands for H±). To further check this conclusion, we will make use of an expansion of 
the cross-sections that is valid in the asymptotic high-mass regime we consider here. The 
cross-sections are simultaneously expanded in m 2 ^ z h /m? Ho and in ^H ,A ,H c VQ/m 2 Ho (except 
maybe for the top quark, the corrections induced by the fermion masses of the SM are really 
negligible in this regime) . The orders of magnitude of these parameters which give the correct 
relic abundance will serve as an a posteriori justification for the use of this expansion, as we 
will see. We separate the various inclusive (co) annihilation cross-sections into A independent 
and A dependent terms as 

^ = o^ + 4 , (28) 

with {i,j = 0, 1, 2, 3} corresponding to {Hq, Aq, H + , H~}. Only the expressions of the dom- 
inant velocity independent terms will be given here. To leading order, for Oq , we have 

11 _ /122 _ j 34 (l + 2<4,)g 4 



411 _ — 

J±n — /In — -ri-i 



USncim 2 



4I3 _ 4 14 _ 423 _ 424 _ s wv 1 rjq\ 



Ho 

2 „4 



647TC? 



rn 



2 
H 
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The A dependent cross-sections can be written in a compact way as 



CT A=^ZI2- . (3°) 



327rm^ o 



with the coefficients A 13 = A 31 given by 



A 00 = A 11 = 2{\* Ho + \\ Q + 2\l c ) 

A 22 =A 33 = 2A 01 = 2(A„ -A Ao ) 2 

A 02 = A 03 = A 12 = A 13 = (X Ho - \ Hc f + (\ Ao - \ Hc f 

A 23 = (X Ho + X Ao ) 2 + AX 2 Hc . (31) 

As we can see, the cross-sections a %3 can only increase when the scalar quartic couplings are 
switched on. The values of the A corresponding to a constant cross-section a 00 (or a 11 ) lie 
on an ellipsoid. For the others a 13 , the ellipsoid is degenerate in a cone (A 02 or A 23 ), or m a 
plane (A 22 ). We can notice in Eq. (|31|) that the coannihilations cross-sections cr^ 1 , <t^ 2 and 
a4 2 are determined by the mass splittings between Hq, Aq and H c , while the annihilation 
cross-sections cj^°, ct^ 1 and cr 23 also depend on the splitting between Hq, Aq, H c and the scale 
111- 

From the positivity of the coemcients A 13 in Eq. (|3ip , we can expect that the relic density 
will decrease when the quartic couplings are turned on. As shown by the instantaneous freeze- 
out approximation, the final relic abundance is actually controlled by the effective thermal 
cross-section Eq. (|16p . where Boltzmann suppression factors e~^ mi ~ m °^ T appear when the 
mass splittings differ from zero. The net result between this thermal damping effect and the 
rise of each cross-section with A turns out to be positive. Therefore, even in the presence of 
non zero scalar couplings, the lower bound rriH > m* for the relic density, with 

m* = 534 ± 25 GeV (3a) , (32) 

remains valid. Above this threshold, the scalar coupling contribution to the cross-section 
has to become progressively dominant over the gauge one in order to obtain the correct relic 
density set by WMAP. It can be seen from Fig. [T] or Eq. (|22p that both contributions become 
equal for mjj ~ y/2m* = 755 GeV. For mu Q > 1.7 TeV, the gauge contribution falls below 
10%. 

For a given mass m# , it is clear that the values of the quartic couplings that are com- 
patible with WMAP are bounded. As we will see, due to Eq. (j3Tj) . they form approximately 
an ellipsoid in the parameter space {Xh , Xa , Xh c }- Mass splittings are also limited, because 

Arriij = mi — rrij ~ (m 2 — m 2 )/2/i2 °c (Aj — Xj) . 

Upper bounds for each |A| and for each mass splitting are shown in Fig. [21 as a function of 
itih - We see that large values of A ~ 0(1) are permitted, the upper bound for each A (or 
a linear combination of them) grows linearly with the DM mass when rriH 3> rn* . We have 
noticed earlier that the coannihilation cross-sections are determined by the mass splittings. 
It turns out that coannihilations involving the charged component are stronger for a given 
mass splitting. This explains why the maximum mass splittings between H c and Hq (or Aq) 
are slightly smaller than the maximum splitting m^ — m jj . Also the maximum rrtn c — mH 
is larger than the maximum m,A — "^h c because we have assumed that Hq is the lightest DM 
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Figure 2: Maximal scalar quartic couplings (left panel) and mass splittings (right panel) as a 
function of mass, imposed by the WMAP bound. Notice that |A# | is maximal for negative 
values of \h - Asymptotic values of the maximum splittings are given in Eq. i33\) . 

particle. As Am^ ~ (m? — m^)/2/i2 ^ (Aj — Xj)vQ/2mH and Aj — Xj oc m# for » m*, 
each mass splitting is bounded by an asymptotic value. Numerically, we find 

\rriA — m-Hol < 16.9 GeV , 
\ m H c — m H \ < 14.6 GeV , 

\rriA — rriH c \ < 13.6 GeV . (33) 

These small splittings serve as an a posteriori justification of the pertinence of the joint ex- 
pansion in rriyy z h // m H anc ^ * n ^h ,a ,h c Vq / we made. They also imply that corrections 
to EWPT observables are negligible in the high mass regime of the IDM (see Eq. ([5])). 

On Fig. [31 two sections in the parameter region allowed by WMAP are shown for three 
values of the DM mass tuhq = 600, 1000, 3000 GeV, they correspond to the two cases itia = 
mH and m// c = {mA +w# )/2. The contours were obtained using the instantaneous freeze- 
out approximation and for a la variation of S^dm h 2 . We have checked that the results agree 
very well with the output of MicrOMEGAs. Ellipsoidal contours are superimposed on these 
figures in red dashed line. They correspond to an expansion in terms of Aniij/T up to 
quadratic terms of the Boltzmann exponential factors in the effective thermal cross-section. 
This ellipsoidal approximation is not accurate when the mass splittings are not negligible 
compared to the temperature around freeze-out, as can be seen for example in the right 
panel of Fig. [3] when m# = 600 GeV. 

We can conclude that the relic abundance required by WMAP can be naturally achieved 
in the large mass regime of the inert doublet model. A viable DM candidate with a mass 
in the multi-TeV range only requires at least one of the scalar quartic couplings to be of 
order 1. In this case, there is no need for any fine-tuning in the parameters of the lagrangian. 
Moreover, the relic density constraint does not put a lower bound on any of the quartic 
couplings (or a linear combination of them). For any value of the DM mass above the limit 
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Figure 3: Contours of X for the WMAP value n DM h 2 = 0.1131 ± 0.0034 for m Ho = 600 
(interior), 1000, 3000 (exterior) GeV, with ijia = mH (left panel) and mu c = (m# + 
m J 4 )/2 (right panel). Dashed curve corresponds to the approximate ellipsoid. 



of Eq. (|25p . one or two of them can even be zero accidentally or if some symmetry is added 
to the model. 

The growth of the values of A with mass as needed by WMAP imposes an upper bound 
on the mass of the DM candidate because of unitarity considerations. If we require all the 
physical quartic couplings Xh , Xa and Xh c in Eq. ([2]) to be smaller than 4ir, we get the 
bound 

m Ho < 58 TeV , (34) 

while for 2ir instead of 4tt we get itih < 30 TeV. This is in agreement with the general 
unitarity bound which holds on any thermal DM relic whose relic density proceed from the 
freeze-out of its annihilation [41]. 

3.2.3 Freeze-out during the unbroken phase of the Standard Model 

We have shown that the inert doublet model in the high mass regime can provide a viable DM 
candidate for a very large range of the DM mass. However, for large masses above ~ 5 TeV, 
the freeze-out will occur before the onset of the electroweak phase transition. In this case, all 
DM components have the same mass mo = /J2 at the epoch of freeze-out, so that deviations 
from the previous analysis are expected. They annihilate into components of the usual Higgs 
doublet Hi in the unbroken phase. The threshold between the unbroken and the broken 
phases occurs at a temperature T c ~ 200 GeV. Although the electroweak phase transition in 
the SM is second order the phase transition, it occurs rather quickly (see e.g. [42] and Refs. 
therein) and for simplicity we will assume a sharp threshold for the freeze-out of our DM 
candidate. This means that the freeze-out will be assumed to be in the broken (unbroken) 
phase for m# < (>)5 TeV. 

In the unbroken phase, the scalar coupling part of the effective annihilation cross-section 
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relevant for the relic density is modified aa^| 



647rmQ 167rmQ ' 

the pure gauge part of the cross-section is still given by Eq. (|22p . in the limit s w — > 0. The 
WMAP constraint determines r\, we get 

r x ~ 2.85 f m ° ) (36) 

Note that r\ = 4tt for mo = 46 TeV and ta = 2n for mo = 22.5 TeV. The DM mass range is 
therefore slightly reduced. 

The maximal values of the scalar quartic couplings and the mass splittings allowed by 
WMAP corresponding to a freeze-out in the unbroken phase are shown in Fig. [U We see 
that |Ah | and A,4 are increased while \jj c is reduced. However, if the stability conditions 
are taken into account, the maximum of |A# | is given by the positive branch A# > which 
is smaller. Also, the maximum mass splitting \mn — m A \ can be slightly higher if the 
freeze-out occurs during the unbroken phase of the SM. We get 

\rriA — mH \ < 17.6 GeV . (37) 

Finally, as it can be seen on Fig. [H the vaccum stability conditions affect significantly the 
maximal values of | m,A — mH | and | rriH c — mu 1 • It is worth emphasizing that the stability 
conditions do not constrain the mass range of the DM candidate. To fulfill them, it suffices 
for A3 to be positive and larger in absolute value than A4 + A5. Therefore, these conditions do 
not put a stringent constraint on the possibility of having a very heavy DM candidate with 
the correct relic abundance. 

3.3 Higher multiplet case 

3.3.1 (Co) Annihilation cross-sections 

In the following we consider only the case of real multiplets. As explained in section 12.21 
for complex multiplets one just must divide the mass obtained with a real multiplet by \pl. 
At tree- level, all the components of a real multiplet have the same mass mo- As the mass 
splittings induced at one-loop are very small (see Eq. ([12]) ), it is a very good approximation to 
consider all states as exactly degenerate. The effective cross-section used for the calculation of 
the relic density is therefore the average of all annihilation and coannihilation cross-sections 
between the odd particles composing the multiplet. The Feynman diagrams of all the relevant 
processes are depicted in the figures of Appendix IB. 21 

For higher multiplets, the only scalar quartic coupling that has an influence on the relic 
density is A3. As for the inert doublet, we can develop the cross-sections in a simultane- 
ous expansion in 'friy/ z h f I m ^ anc ^ m ^3 v o/ m o- At leading order, the dominant velocity 
independent terms of the effective annihilation cross-section a^v = (Tq^v + cr^v ar^| 

4 „, = („»-!)(,,' -3) I md A W_I_a, . ( 38 ) 

n 128vr ml x n 16tt mg 



3 We assume that all the components of the Higgs doublet have masses much smaller than mo- 

4 Like in the doublet case, the expression of A n ^ in Eq. (|38f) differs from the result of Ref. [25] by a factor 



1/2. 
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Figure 4: Maximal values of scalar quartic couplings (left panel) and mass splittings (right 
panel) as a function of the DM mass, constrained by WMAP, without (dashed lines) and 
with (thin solid lines) the vacuum stability conditions Eq. ^ included. We assume a Higgs 
mass rrih = 120 GeV, and a sharp threshold between the freeze-out in the broken and in the 
unbroken phases of the SM at a mass m# = 5 TeV. 



They drop as m,Q 2 , as expected from unitarity considerations. In the pure gauge limit 
A3 = 0, the odd DM components annihilate almost exclusively into gauge bosons. Coan- 
nihilations into fermion final states are p-wave suppressed. When A3 7^ 0, new channels of 
(co) annihilation are opened, through a Higgs particle, or into Higgs particles. Expressions for 
the velocity dependent terms in the cross-section will not be given, as they are subdominant. 
We have checked numerically that they lead to a maximal correction smaller than about 5%. 
They have been taken into account in a numerical evaluation of the relic density with the 
instantaneous freeze-out approximation. 

Notice that for high multiplets (n > 2), the high mass regime we consider (moM > m w) 
is the only possibility for a successful DM phenomenology. Below mw-i given the collider 
bounds on the charged multiplet component and given the small neutral-charged component 
mass splittings, coannihilation cross sections would be far too large to account for WMAP 
DM abundance. 

3.3.2 Relic density 

The relic density in both real and complex models with n = 3, 5, 7 has been computed using 
MicrOMEGAs [39], and compared to the result of the instantaneous freeze-out approximation. 
The agreement between the two approaches is better than 4.5 %. For a real multiplet of a 
given dimension n, the relic abundance SIdm^- 2 depends only on the two free parameters of 
the model, the mass of the DM candidate mo and the coupling A3. Therefore the WMAP 
constraint on the relic density determines A3 as a function of mo, or vice- versa. The values 
of mo corresponding to A3 = 0, 2ir, 4-7T are given in Table [TJ We find threshold masses (i.e. 
for A3 = 0) that are systematically smaller than the values obtained in Ref. [25] by ~ 10%. 
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Models 


A 3 = 


A 3 = 2?r 


A 3 = 47T 


A 3 = (SE) 


A 3 = 4?r (SE) 


Real Triplet 


1.826 ±0.028 


11.1 


21.9 


2.3 


28.1 


Real Quintuplet 


4.642 ±0.072 


9.6 


17.4 


9.4 


35.7 


Real Septuplet 


7.935 ±0.12 


10.6 


16.1 


22.4 


46.3 



Table 1: Threshold masses (in TeV) without or with Sommerfeld effect (SE) for higher mul- 
tiplet models, as determined by the WMAP constraint, the errors quoted correspond to a la 
variation of the relic density. The large mass range of the DM candidate is shown by the 
indicative values for A3 = 2ir and 4ir. 

The values of the parameters niQ and A3 that are in agreement with the WMAP constraint 
at 1 a level are shown on Fig. [5] for all real multiplet models. Similarly to the doublet case, 
the following mass-coupling relations hold: (mo — m^) oc A| for tuq close to the threshold 
value mj, while A3 oc mo for mo S> mj. This behaviour is easily recovered from the analytic 
expression of the effective cross-section, Eq. (|38p . because the WMAP constraint fixes its 
value. More precisely, we see that for mo 3> mj, the DM mass scales like mo ~ which 
explains the different slopes of the linear part of the function mo (A3) (see Fig. [5]). 

An upper bound on moM can be obtained by demanding that the theory stays perturba- 
tive. Values of mo corresponding to A3 = 2ir and A3 = 4tt are given in Table [TJ For higher 
multiplet models, the allowed DM mass is in the mufti- TeV range, even when the Sommerfeld 
corrections are not included. For all candidates with a mass higher than around 5 TeV, the 
freeze-out will occur in the unbroken phase of the SM. Unlike the doublet case, the expres- 
sions for the effective cross-sections given by Eqs. (|38p for the broken phase remain valid in 
the unbroken one, although the detailed (co) annihilation processes are different. Therefore, 
the behaviour of mo as a function of A3 given in Fig. [5] is still valid. For higher multiplets, 
the so-called Sommerfeld effect plays however a significant role. 

3.3.3 Sommerfeld effect 

At small relative velocity, the interaction between two particles becomes long range if the 
mass of the particle exchanged in the interaction is much smaller than the two interacting 
particles masses. This leads to a non perturbative enhancement of the annihilation cross- 
sections of very heavy DM candidates, known as the Sommerfeld effect [43]. For an abelian 
vector interaction, the two parameters which determine the strength of the enhancement are 
a/e = m DM /(rnv/a) and a/ (3, where my is the mass of the vector particle, a is the coupling 
constant and (5 is the DM velocity. For non abelian interactions, the general trend of the 
enhancement is controlled by the same parameters than in the abelian case, but is complicated 
by the possibility of resonances [44]. It has been shown in Ref. [28] that Sommerfeld effect 
corrections affect both the relic density calculations and the present day annihilation cross- 
sections relevant for indirect detection signals. 

A full treatment of the Sommerfeld corrections is beyond the scope of the present pa- 
per. In the case of heavy scalar dark matter candidates, corrections to both weak and 
scalar interactions are a priori expected. However, given the fact that the scalar interactions 
give contributions either pointlike or, from the exchange of a Higgs boson, suppressed by a 
v o/ m DM factor, in what follows, we will only take into account the Sommerfeld corrections 
due to gauge bosons, which has been calculated in details in Ref. [28]. As their full computa- 
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Figure 5: Evolution of the mass of the dark matter candidate as a function of the coupling A3 
for all the higher multiplet models of phenomenological interest, as constrained by WMAP, 
without (solid lines) or with (dashed lines) Sommerfeld effect. The curves correspond, from 
top to bottom at A3 = 0, to the real septuplet, the real quintuplet and the real triplet. The 
shaded area on the left is excluded by the vacuum stability constraint (for = 120 GeV and 
\™ ax = 4tt). 



tion shows, away from resonances, the enhancement factor at the time of freeze-out is almost 
constant. Moreover, this factor strongly depends on the multiplet dimension. For the dou- 
blet, it is negligible (a few percent), but it increases rapidly for higher multiplets. The first 
result can be understood by the fact that thermal contributions to the gauge boson masses 
were included. If my oc gT, the parameter a/e is roughly constant at the time of freeze- 
out T ~ tudm/xf- The second result can be understood by the fact that (co) annihilation 
cross-sections grow roughly as n 2 with the dimension of the multiplet (see Eq. (|55p for the 
annihilation of Ao). Effectively, the coupling constant is therefore enhanced by the size of 
the multiplet. Numerically, we find that n 2 — 1 = 8, 24, 48 for n = 3, 5, 7 scales well as the 
enhancement factor ~ 1.6, 4.1, 8.0. 

In Fig. [5l we show how the Sommerfeld effect modifies the relation between mo and A3 
for all higher multiplets. Following the line of reasoning of the last paragraph, a simple 
approximation has been used. A constant cross-section enhancement factor is taken for each 
multiplet. It is chosen so as to reproduce the threshold masses given in Ref. [28] (see Table [Q) 
when Sommerfeld corrections are taken into account. Resonances as well as the possible 
enhancement factor from scalar couplings have been neglected in this simple treatment. Such 
effects would push the DM mass to even higher values for a given value of A3. 

4 Direct detection 

The neutral scalar field of an SU(2)l multiplet has vector-like interactions with the Z boson 
if its hypercharge Y 7^ 0. This leads to an elastic spin independent cross-section between the 
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DM candidate and the nucleon that is 2 to 3 orders of magnitude above current limits [45,46]. 
To evade this constraint, either the multiplet has Y = 0, or a mass splitting (at least of the 
order of ~ 100 keV) is induced by some mechanism between the real and the imaginary 
components of the neutral field. In the latter case, the DM - nucleon interaction through 
the Z boson is kinematically forbidden, or leads to tiny inelastic collisions. The doublet 
case is special in this respect, because the most general renormalisable potential Eq. ([I]) 
automatically allows for such mass splitting. This direct detection constraint explains why 
we only considered the inert doublet model and higher multiplets with n = 3,5,7 in this 
work. 

For the models studied, the only tree-level interaction between the DM candidate and 
the nucleon proceeds through the exchange of a Higgs scalar, which gives rise to elastic spin 
independent collisions. If is the coupling between the DM and the Higgs particle, the 
DM-nucleon cross-section is then given by 

^MN*/&f(-^Y . (39) 

where /jv — 0.3 is the nucleonic form factor determined experimentally (see [47] and also, 
for a range for /jy, see [13] and references therein.). For the doublet case, = A# while 
Ah = A3/2 for higher multiplets. 

Elastic scatterings through the exchange of SU(2)l x U(1)y gauge bosons is also possible, 
but only at loop level. Its contribution to the elastic cross-section has been computed in 
Ref. [25], 

4n = fUn 2 ~ If^f- Mr + ^Y (40) 
64m^ \mfo m hJ 

The value of this pure gauge part increases rapidly with the dimension n of the multiplet, 
and does not depend on the DM mass for a given multiplet. 



4.1 Doublet 

The mass splitting between #0 an d Aq is controlled by A5. As we have seen, the relic density 
constraint does not put a lower bound on the absolute value of this parameter. Tiny mass 
splittings are allowed, and are stable against radiative corrections. This opens the interesting 
possibility of explaining the DAMA annual modulation data while still being compatible with 
other experimental bounds through inelastic collisions with the nucleon Hq n — > Aq n [48]. 
It has been shown that the mass splitting \rriA — ?Ti# | should be roughly in the range 
[50. ..150] KeV to realize this scenario (sec e.g. [49] and Refs. therein). For a DM mass 
between 535 GeV and 10 TeV, this corresponds to a tiny value of A5, 

0.9 • 10~ 6 < A 5 < 5 • 10~ 5 , (41) 

a range that also has important possible implications on leptogenesis and neutrino masses 
(see Section [6] below). A precise determination of the parameter range that fits all direct 
detection data strongly depends on assumptions on the velocity distribution of the dark 
matter particles in the Earth neighborhood, and will not be presented here [48]. In any case, 
the inelastic collisions rapidly become inoperant when the mass splitting is increased above 
1 MeV. 
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Figure 6: Elastic cross-section on nucleon for the inert doublet (left panel) and for higher 
multiplets (right panel), compared to experimental limits (CDMS Ge result from 2008 [45], 
Zeplin III final result (2008) [46]) and projected sensitivities at future experiments (Super- 
CDMS and Xenon IT) [50]. We have assumed = 120 GeV, a standard Maxwellian 
DM halo with a local density po = 0.3 GeV/cm 3 . For the left panel the shaded area gives 
the allowed range of values. Its lower limit corresponds to the pure gauge interaction cross 
section putting all quartic interactions to 0. The upper limit on the elastic cross-section is 
given by the solid (dashed) blue line when vaccuum stability conditions are (not) taken into 
account. For the right panel, solid (dashed) curves correspond to the cross-section prediction 
without (with) Sommerfeld effects. 



In the case where the inelastic collisions through a Z boson can be neglected, we are left 
with the elastic cross-sections of Eqs. (|39M40p . The maximum value of Xh allowed by the 
relic density constraint grows linearly with mjy , this translates into an absolute upper bound 
on the elastic cross-section. Numerically, we find 

<t HoN < 9.4- 1(T 9 pb . (42) 

While on one hand an upper bound is derived from the WMAP constraint, on the other 
hand, the pure gauge cross-section Eq. (|40p sets a lower bound around 10~ 10 pb. As a result, 
the direct detection rate can vary by two orders of magnitude, for a given DM density and 
velocity distribution around the earth neighborhood. In Fig. EJ we show the range of the 
elastic cross-section as a function of mass, compared to current experimental bounds and 
future experiments sensitivity. The DM-nucleon cross-section has been calculated with m/i = 
120 GeV. The limits and projections assume a standard value for the local DM density po = 
0.3 GeV/cm 3 , and a maxwellian velocity distribution with the characteristic halo velocity 
v haio = 220 km/s. Extending the expected XENON reach to 10 TeV, we see that a large 
portion of the parameter space of the IDM in the high mass regime will be probed by future 
direct detection experiments with a 1 Ton x year sensitivity. 
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4.2 Higher multiplet 

In the case of higher multiplets, no mass splitting between the neutral components of a 
complex multiplet can be generated with the scalar potential of Eq. 0. As a result, only 
multiplets with vanishing hyper charge Y = are viable. Therefore, we are led to consider 
only real multiplets of dimension n = 3, 5, 7. For a given DM mass, the only free parameter is 
A3, which is determined by the relic density constraint (see Fig. [5]). Higher multiplet models 
are therefore particularly predictive. 

As for the doublet case, the spin independent elastic scattering cross-section has two 
parts. The pure gauge part given by Eq. (|4U|) yields 



The scalar quartic coupling part increases with ttt-dm, with the following upper bounds, 



The total cross-section is represented in Fig. [6l It is well above the sensitivity limit of the 
future XENON IT experiment but below current limits. Contrarily to the doublet case, 
the gauge contribution to the elastic cross-section is always substantial or even dominant 
compared to the scalar one. When Sommerfeld corrections to the relic density calculation 
are taken into account, the relative contribution of the scalar interactions to the elastic cross- 
section is even smaller (see Fig. [6]) . In all cases for a given value of m dm the direct detection 
cross section is predicted to one value. 

5 Indirect detection 

The annihilation of DM can produce several types of signals useful for indirect detection 
searches. We will examine photons and neutrinos, which give a directional signal, and also 
charged antimatter cosmic rays for which such a directional information is lost after diffusive 
processes. Heavy scalar DM particles mainly annihilate into ZZ, W + W~ and hh, which 
subsequently produce the desired signal in cascade decays. Therefore, the annihilation of 
heavy scalar candidates generally produces soft spectrums. 

In what follows, we will only derive predictions for these soft spectrums. The monochro- 
matic signal from direct annihilation into photons at one loop for example will not be consid- 
ered here, as this cross-section is strongly affected by non perturbative effects when it becomes 
non negligible [51]. Generally speaking, a detailed analysis of Sommerfeld enhancements and 
resonance effects is beyond the scope of the present paper, so that all predictions will be 
made with an enhancement factor EF = 1. It is however worth noticing that heavy scalar 
DM models considered here are viable for a wide range of mass when scalar interactions are 
taken into account. As a result, the phenomenon of resonances, which occurs for particular 
values of the DM mass, can always be achieved by tuning the DM mass to one of these values. 
Such a possibility does not appear for fermionic minimal dark matter candidates, where the 
DM mass is determined by the relic density constraint [25]. 

We will now give a brief description of the flux calculation for each type of indirect signal. 
Then, the predictions for the inert doublet and the higher multiplet models will be presented. 



cr° MN = (0.86, 7.76, 31.03) x 10" 9 pb for n = 3, 5, 7 . 



(43) 



°dmn < ( 3 - 34 > 5 - 28 > 6 - 15 ) x 10 P b for n = 3, 5, 7 . 



(44) 
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5.1 7 an d v signals 

The galactic center (GC), where the DM halo surrounding the galactic disk is believed to 
be the most concentrated, is the most promising region to probe for DM annihilations. If 
the halo is very cuspy, as simulations [52-54] and some dynamical mechanisms [55] suggest, 
observation of photons or neutrinos from the GC can provide a clean signature of the presence 
of DM. However, if the halo has a rather flat profile, as indicated by direct kinematical obser- 
vations [56] , the signal might be difficult to disentangle from the astrophysical background. 
The total flux of 7 or v in a solid angle Af2 around the galactic center is simply calculated 

as 

^ (A0 ) = JqL at x An f° Ro J(Afl) , (45) 

DM 4?r 

where 

is the average number of 7 or u per annihilation with an energy between the experimental 
thresholds E m i n and E max , po is the local DM density, and 

dc 1 r 

J(ASl) = 2 / P 2 dldn, (47) 

is a dimensionless astrophysical factor which encodes all the uncertainties about the distribu- 
tion of DM in the galactic halo. The quantity BF is the so-called (astrophysical) boost factor, 
an enhancement factor due to the dumpiness of DM in galactic halos. It should however be 
stressed that the boost factor is dependent upon the observation direction. For the direction 
of the galactic center, BF is negligible if the halo is very cuspy. For a flat profile like the 
isothermal one, BF is also limited because the concentration of subhalos is comparable to 
that of the galactic halo [57]. On the particle physics side, the Sommerfeld effect can provide 
non negligible enhancements. The predictions made hereafter do not include any boost. 

For the sake of completeness, let us finally mention that in the case of neutrinos, an 
interesting possibility is to search for an signal from the core of the Sun or of the Earth, 
emanating from annihilations of DM particles captured by the celestial body. However, in 
the high mass regime, the capture rate of DM particles, which scales ELS ^^JJfyJ IS 

too small to 

lead to observable signals. Also, in the case of the Earth, it cannot be enhanced by resonance 
effects like for lighter candidates with a mass around 50 GeV [32]. 

5.2 Charged antimatter cosmic ray signals 

As antimatter cosmic rays (CR) are quite rare in the galaxy, they are also promising mes- 
sengers to probe for exotic physics, like DM annihilations. The recent publication of the 
positron fraction observed by PAMELA [58], together with the excess seen by the ATIC ex- 
periment [59] have triggered a lot of activity in this research field, as they point to a positron 
excess between 10 and 800 GeV [60]. The explanation of both excesses by a DM scenario 
would lead to a candidate with rather unusual properties [61], and is therefore not favored. 
In this paper, typical flux spectrums for both positrons and antiprotons are presented, but 
no attempt will be made to fit the PAMELA or the ATIC excess. 



21 



The diffusive random walk of charged charged cosmic rays through the galaxy can be 
described by the following general (steady state) propagation equation [62]: 



V 



K{E,x)VM cr -V com M a 



+ 



d_ 
~dE 



b(E)M cl . + K EE -^M a 



+ T(E)M CT + Q = (48) 



where jV cr = dn cr (E)/dE is the number density per unit of energy, K(E,x) and K EE are 
coefficients that encode the diffusion by galactic magnetic fields in real and momentum space, 
V conv is the velocity of the galactic convective wind, b(E) is the rate of energy loss, T(E) 
accounts for spallation processes (destruction of CR due to collisions with the interstellar 
medium), and Q is the source term. For positrons, the dominant processes are the energy 
loss and spatial diffusions. For antiprotons, the dominant processes beside diffusion are 
spallation and convection. For a detailed discussion of the propagation model, we refer the 
reader to Ref. [63-66]. 

The flux of a cosmic ray species cr at the earth location is obtained by convoluting the 
Green function of the propagation equation with the source term Q, given by 



Q = BF 



(av)p 2 



2m 



DM 



(49) 



It is worth emphasizing that the astrophysical boost factor in this equation is in general 
energy dependent. The imprint of the dumpiness of the DM halo on the CR spectrum indeed 
depends on the typical diffusion length, which itself depends on the injection energy. It has 
been shown [67, 68] that only high energy positrons and low energy antiprotons (compared 
to the injection energy) can be sensitive to a variation of the local DM density. For the rest 
of the spectrum, the astrophysical BF never exceeds one order of magnitude. 

In this work, the propagation has been carried out by the code DarkSUSY [69]. The 
(default) propagation models implemented in this package correspond to a simplified version 
of Eq. (|48j) . where only the most relevant processes specific to each CR species are included. 
Finally the effect of the solar wind on charged particles is taken into account by applying the 
force-field approximation, which results in a shift in energy between the interstellar spectrum 
(IS) and the one at the top of the atmosphere (©), 



Eis = + \Ze\4> , 
and a depletion of the flux at low energies (below ~ 10 GeV) 

dE® p 2 IS dE IS 



(50) 



(51) 



where and pjs are the momenta at the Earth and at the heliospheric boundary, Ze is 
the charge of the CR particle and cf> is the solar modulation electric potential which we took 
equal to 600 MV (see e.g. [70] and reference therein.). 
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Point m Ha (TcV) av (pb) BR(ZZ) (%) BR(W+W-) (%) BR(fefe) (%) BR(tt) (%) 



I 


1.0 


1.47 


21.0 


24.8 


49.7 


4.5 


II 


1.0 


1.56 


77.7 


22.3 








III 


1.0 


1.76 


17.7 


82.3 








IV 


10.0 


1.28 


0.2 


0.3 


99.4 


0.1 


V 


10.0 


1.35 


99.75 


0.25 








VI 


10.0 


1.64 


0.2 


99.8 








VII 


10.0 


3.0 


25.0 


50.0 


25.0 






Table 2: Annihilation cross-section and branching ratios of the inert doublet DM candidate 
in the high mass regime. These benchmark points give the correct relic abundance as required 
by WMAP. 

5.3 Doublet 

The annihilation cross-sections of Hq into ZZ, W + W~ , and hh pairs are given by 

a{H H ^ZZ)v = — \ 2 + A ° 



1287rcim 2 Ho 16vrm^ o 
a(H H ^W + W-)v = —4-^ + 



64irci<m 2 Ho 8vrm^ o 
A 2 

a(H H ^hh)v = H \ . (52) 
\bTTm Ho 

To estimate the various indirect detection signals, we consider seven benchmark points that 
lead to the relic density required by WMAP Each of the first six points maximizes the 
branching ratio in one of the three possible dominant annihilation channels of Hq, for two 
values of the DM mass m# = 1 TeV and m# = 10 TeV. These benchmark points serve to 
evaluate the spread due to parameters on the particle physics side. However, the point IV 
does not satisfy the stability conditions of the potential. If these were taken into account 
(point VII), the maximum branching ratio into hh would be around 25%. For this last 
point, the annihilation cross-section is significantly higher. It corresponds to the degenerate 
limit situation ~ mA — ht-Hc where coannihilations during the freeze-out epoch are 
the strongest. Table [2] gives the annihilation cross-section and the branching ratios for the 
benchmark points. 

The numbers of photons and of neutrinos from the galactic center can be calculated with 
Eq. (I45p . We will estimate the flux from a solid angle AQ = 10~ 3 corresponding to a cone 
with an aperture of 2°. For photons, a typical experiment like FERMI-LAT (former known 
GLAST) has an angular resolution 5Q ~ 10~ 5 and energy thresholds 1 < E~ < 300 GeV. 
For mH = 1 TeV, the number of photons per annihilation 25 < A 7 < 40 is slightly higher 
for annihilations into hh. For m# = 10 TeV, 30 < iV 7 < 100 is highest for annihilations 
into If we consider a cuspy halo profile like the NFW (Navarro- Frenk- White) one, 

J ~ 1.3 • 10 3 , which gives 

$ 7 (Aft = HT 3 ) ~ Oil) x 2.3- 10~ 10 (^^X 2 [phcm^s- 1 ] , (53) 

\ 1 TeV J 

which has to be compared to the FERMI-LAT sensitivity at about 10 -10 phcm~ 2 s~ 1 for 
point sources [71,72]. For a flatter profile like the isothermal one (J ~ 25), the signal would 
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Figure 7: Positron (left panel) and Antiproton (right panel) flux resulting from H$ anni- 
hilation for the benchmark models I-VI of Tabled The data are taken from [59, 7^-78] 
and from [79-83]. These figures were obtained with a standard NFW halo profile, a solar 
modulation potential (f> = 600 MV and no boost factor. 



be broadly distributed over the bulge region. Even if the total flux lies above the sensitivity 
for a diffuse flux at about 10~ 10 phcm" 2 s" 1 sr" 1 [72], it would be difficult to disentangle the 
DM signal from the astrophysical background. 

For neutrinos, we consider the experiment Antares and the forthcoming extension KM3net, 
for which the galactic center is visible. The energy threshold is E u > 100 GeV, and the typ- 
ical angular resolution for 1 < E v < 10 TeV is 5Q ~ 10 -4 [73]. The number of neutrinos N v 
is suppressed in the case of annihilations into hh and their energy spectrum is softer. For 
annihilations into ZZ or W + W~, N u ~ 10 (1.0) for m Ho = 10 TeV and E™ in = 100 GeV 
(1 TeV). Using again a NFW profile, we get a flux 

mfcev) t^ 111 " 2 ^ 1 ] . ( 54 ) 

while the point source sensitivity of KM3net for 1 year will be ~ 3 • 10~ 10 z^cm _2 s _1 (3 • 
lO" 11 i/cm^s" 1 ) for E'f m = 100 GeV (1 TeV) [73]. Therefore, the detection of neutrinos 
from annihilations of our DM candidate in the galactic center looks very difficult. 

Positron and antiproton fluxes for the benchmark models of Table [2] are shown in Fig. [7] 
and we see that they lie well bellow the data. For m# = 1 TeV, a boost factor of ~ 3 orders 
of magnitude would be needed to reach the range of the observed positron flux. In our result, 
the boosted antiproton flux would however still be below the background. In the framework 
of the inert doublet model, such a boost factor is possible only if the DM mass is very close 
to a resonance. The detailed analysis of this possibility is beyond the scope of this paper. 
Going to higher masses decreases the number density of DM particles in the halo. The CR 
fluxes produced by the annihilation of DM are therefore even more suppressed for a DM mass 
larger than 1 TeV. For rriH = 10 TeV, the positron flux is at least 4 orders of magnitude 
below the background signal. Predictions for the positron fraction in the IDM, as well as a 
comparison with the excesses seen by PAMELA and ATIC can be found in Ref. [33]. 
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Model 


av (pb) 


BR(ZZ) (%) 


BR(W+W~ 


) (%) BR(ftft) (%) 


Real Triplet 


2.47 


24.4 


51.1 


24.4 


Real Quintuplet 


3.81 


21.7 


56.5 


21.7 


Real Septuplet 


4.19 


13.1 


73.7 


13.1 



Table 3: Annihilation cross-section and branching ratios for higher multiplet candidates with 
a mass mo = 10 TeV and a relic density set by WMAP. Sommerfeld corrections are omitted. 



Finally, we can notice that the CR spectrums are significantly softer in the case of an 
annihilation into a pair of Higgs particle. The Yukawa coupling of h to fermions is directly 
proportional to the fermion mass, and becomes small compared to gauge couplings for light 
fermions. Higgs particles decay mainly into bb pairs, leading to multiple hadronization cas- 
cades and jets. The enhancement of the antiproton flux at low energy is particularly clear (1 
order of magnitude!). 



5.4 Higher multiplet 

In the case of higher multiplet models, the annihilation cross-sections into ZZ, W + W~ , and 
hh pairs are given by 



ct(A A -> ZZ)v 



647rmQ 



(n 2 -l)V , A| 



ct(A A -»• W + W~)v = v ncr , ' a + 



2567r?T7,Q 327T7T7.0 

ff(AoA -WO« = 77^4 • (55) 

With a large gauge contribution, the W + W~ annihilation channel is always dominant. Omit- 
ting the Sommerfeld enhancement, we give in Table the typical annihilation cross-section 
and the branching ratios for a candidate with a mass niQ = 10 TeV, and the required relic 
density for n = 3, 5, 7. Although the annihilation cross-section increases with the multiplet 
dimension, the conclusions obtained in the doublet case for the detectability of the various 
indirect detection signals still apply. Gamma ray telescopes offer the most promising search 
and a better sensitivity than the neutrino detectors. Again, the production of charged cosmic 
rays is well below background unless a huge boost factor is applied. 



6 Neutrino masses, leptogenesis and DM at a low scale in the 
doublet case 

If, to account for the neutrino masses, one adds right-handed neutrinos iVj to the inert 
doublet model there are 2 possibilities: either these TV's are even under Z2 or they are odd. 
One interesting consequence of the results obtained above for the doublet case, which arises 
for the latter possibility, is that they allow successful generation of both neutrino masses and 
baryogenesis via leptogenesis in a way where DM plays an important role. Moreover all three 
phenomena can be induced at a scale as low as TeV, even for a hierarchical spectrum of iV's. 
This has to be compared with the lower bound which exists on the mass of the right-handed 
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Figure 8: One-loop neutrino mass diagram. 



neutrinos, m^v > 6 • 10 s GeV [84-86], for a hierarchical spectrum of iV's in the usual type-I 
seesaw model where only right-handed neutrinos are added to the SM. To our knowledge, 
the mechanism we propose in this section is the most simple and minimal way to induce all 
three phenomena at such a low scale in a related way. 

The crucial point at the origin of the fact that leptogenesis can be generated at such a 
low scale is that, if the iV's are odd under the Zi symmetry, Yukawa coupling involving the 
iV's and the Higgs doublet are forbidden. The most general lagrangian one can write is [18] 

C = C IDM + WipNi - NiY Nij H\Lj - X -m Ni NM (56) 

with H2 = ^H^, i.e. only the Yukawa couplings with the inert doublet are allowed. As a 
result neutrino masses cannot be generated at tree-level in the usual way but only at one 
loop through two DM inert doublets, Fig. [HJ which for m^i » m H o A o Hc gives (at lowest 
order in As^/m^ ) 

,2 



Mv^^flog^ + ll. (57) 



m N k 



With respect to the standard tree-level seesaw model, which gives (m u )ij = — ^ Xe^EM ; this 
" radiative seesaw" mechanism leads consequently to an extra suppression of the neutrino mass 
by a factor ^^-[log(m^ / 'm% k ) + 1] for each Nj contribution. 

As for leptogenesis in this framework, it proceeds from the N — ► LH2 , LH2 , that is to say 
in the same way as in the usual type-I seesaw model, replacing all ordinary Yukawa couplings 
to a Higgs doublet by the inert doublet Yukawa couplings of Eq. ([56]) . Fig. [9]. For the lightest 
right-handed neutrino N\ and mjv 2 3 » this gives the CP-asymmetry 

^ 3 m^EiMO^w^tj) 2 ] fKQ s 

It is the extra suppression above of the neutrino masses versus absence of any extra suppres- 
sion of the CP-asymmetry which allows to lower the scale of leptogenesis, as we will now 
show in details by deriving the various relevant bounds: 

1. Leptogenesis and neutrino mass bounds on £jv 15 rn^ 1 and A5. In full generality 
in the type-I seesaw model, ejvj is bounded by the size of neutrino masses [86] 

3 mjvj \Im[Trm^ml' 3 ]\ 3 

RVi = o 2 = ^ q 2-{ m v 3 ~ m vi) ( 59 

8tt Vq mi 8-7T Vq 
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Figure 9: One-loop diagrams contributing to the asymmetry from the iVj decay, involving 
the DM inert doublet. 



with mj, the contribution of Ni to the neutrino mass matrix, m U3 (m Ul ) the mass of the 
heaviest (lightest) neutrino, and mi = J2i I YNli\ 2 Vo/2mN 1 . This leads to the lower bound 
mjviJ>6 • 10 8 GeV, imposing that the baryon asymmetry produced ub/s = — Hn^/s = 

~^^r'% £N i r l = ~ L35 ' w ^ 3£ N 1 r] is at least equal to the WMAP value, n B /s = 9 • 1CT 11 
(assuming a maximal efficiency rj = 1), with g* = 108.75 the number of active degrees of 
freedom at the time of the creation of the asymmetry. For the radiative seesaw case, taking 
for simplicity all [log (m^/m^- ) + 1] factors equal to unity, one gets an extra 8-7r 2 /A5 en- 
hancement factor of the CP-asymmetry and hence a decrease of the mjVj lower bound by the 
same amount H 



87r 2 3 m/v, \Im\Trnh) ml' 3 ]\ . 8tt 2 3 m N , 
eJVj = -r-g 2 = ~ "A - i — m ^ " m »i) 60 

A5 87T V z mi A5 87T v z 



m N,>p2 6 - 108 GeV ( 61 ) 

The change of number of effective degrees of freedom due to the extra active inert compo- 
nent^) is of little importance and can be neglected. 

Consequently for successful leptogenesis the lower bound can be lowered down to any 
value mjvi if 

A 5 < 1.5-10 -4 -(T7WlTeV) <->• rn Ao -m Ho < 9MeV-(500GeV/m DM )-(m A r 1 /lTeV) (62) 

2. Bound on A5 from DM constraints in the low mass regime and corresponding 
lower bound on m^i for successful leptogenesis. In the low mass DM regime, where 
771dm < TTHyjj to avoid too fast Hq-Aq coannihilation leading to too low relic density, it is 
necessary that the mass splitting is large enough, mjj — m,A > 7 GeV [19], or equivalently 
A5 > 1.6 • 10~ 2 . This bound is incompatible with successful leptogenesis, i.e. Eq. ([62]) . unless 

m Low regime > nQ (g3) 

which is interestingly low but not enough to be reachable in a not too long term at colliders. 



5 Notice that if one doesn't take the approximation that all logarithmic factors are unity one can get even 
a much more relaxed bound, taking for instance [log(m| f /m^ ) + 1] ~ 0. We will not consider this peculiar 
possibility preferring here to stay generic. 

6 Leptogenesis in this case has been considered in Ref. [87] for large right-handed neutrino masses around 
10 9 GeV. 
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3. Lower bound on m^ L for successful leptogenesis in the high mass regime. In the 

high mass regime, as we have seen in Section. l3,2.2p and unlike in the low mass regime, the relic 
density constraint doesn't lead to any lower or upper bound on A5 (apart from a perturbativity 
bound). As explained above the only lower bound comes from direct detection constraint, 
mA — mHo ^ 100 keV, which is well below the bound of Eq. (|62j) (unless mjvi ^20 GeV, 
that is to say well below the sphaleron decoupling temperature anyway). The only relevant 
constraint in this case comes therefore from successful conversion of the lepton asymmetry 
into a baryon asymmetry by the sphalerons. Assuming a sphaleron decoupling scale of order 
150 GeV [42] and taking into account that the creation of the lepton asymmetry cannot be 
instantaneous, one gets the constraint: 

m Hi9hregime> 1TeV (g4) 

4. Bounds on the iVj Yukawa couplings. All bounds above are obtained assuming that 
there is no efficiency suppressions of the lepton asymmetry produced. This requires that the 
decay width of N\ satisfies the out of equilibrium condition: 



1 Uir 3 a T 2 
= —\Y Nlj \ 2 m Nl < H(T = m Nl ) ~ \ —f- — (65) 



8tt' - 1 y V 45 M Planch 

which gives 



T=m N 



\Y Nlj \ l < 4- 10" 14 • (mjvJITeV) (66) 

Notice that the smallness of these couplings doesn't induce any suppression of the CP- 
asymmetry because these couplings essentially cancel in it (up to phases), see Eq. (|58ft . 
Eq. (|66p implies that Ni gives neutrino mass contributions smaller than the atmospheric 
or solar mass splittings. These splittings must therefore be dominated by the contribution 
of -^2,3 and the neutrino mass spectrum is necessarily hierarchical. The heavier A^2,3 states 
must have necessarily larger Yukawa couplings. Eq. (|62p . together with the lower bound 



m„ 3 > y Am^ ta = 0.06 eV and Eq. flSTJ) imply a lower bound on the N2 or N3 Yukawa 
coupling 

Y Nji > 1 • 10- 3 • 2- ) (j = 2 or 3) (67) 

for at least one lepton flavor i value. This bound can also be obtained directly from Eq. (|58p 
imposing that the CP-asymmetry is large enough to induce the observed baryon asymmetry. 

Comparing Eq. (166p with Eq. (167p . one therefore concludes that to induce leptogenesis 
at a low scale in a non-resonant way [85,88,89] one needs a hierarchical structure of Yukawa 
coupling (similar to the one of the charged leptons). 

Eq. (|67p gives Yukawa coupling values much larger than in the usual type-I seesaw model 
where in order to give neutrino masses ~ \J Am^ tm one needs Yukawa couplings ~ 0(1O~ 6 ) 
for tojv ~ 1 TeV (unless cancellations occur between the Yukawa couplings in the neutrino 
masses) . 

5. AL = 2 washout constraint. Finally in order to have an efficiency of order unity as 
assumed above it is also necessary that there is no washout from scattering processes. The 
most dangerous are the ^2,3 mediated AL = 2 ones, due to the fact that the A^2,3 Yukawa 
couplings must be fairly large. However for values of Yat 2 3 of the order of the bounds in 
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Eq. (|67|) . it can be checked from the Boltzmann equation [85,88,89] that this effect is moder- 
ate or even negligible (even without needing to play with flavor effects which could be invoked 
otherwise to suppress further this AL = 2 washout effect). 

Taking into account all constraints above, as a numerical exampld^, for A5 = 10 -4 , = 
2 TeV, m N2 = 6 TeV, Y Ni ~ few 10~ 8 , max,(Y jv 2 i) 4 • 10~ 3 and m Ho above ~ 510 GeV 
and sizably below mjVj we get m Vz = \J Am^ tm , the WMAP value n^/s ~ 9 • 10~ n (with no 
sizeable suppression of the efficiency) and a dark matter relic abundance which can be easily 
consistent with the WMAP range above (i.e. in agreement with the results of Section. I3.2.2j> . 
An interesting property of this framework is that it involves Yukawa coupling of N2 and/or 
N3 much larger than in the usual seesaw model. Unlike in the latter case, it is therefore 
conceivable in the not too long term to produce the right-handed neutrinos of the present 
extended IDM at colliders. This is an interesting phenomenological possibility to test in 
addition to the nature of DM, the related origin of neutrino masses and baryogenesis via 
leptogenesisJl 



7 Summary 

Properties of scalar DM candidates with SU(2)l quantum numbers are driven by their known 
gauge interactions and by their scalar quartic interactions. If the quartic couplings are not 
much smaller than the gauge interactions their effects cannot be neglected. This leads, for 
each model, to a range of values of DM masses which can reproduce the observed DM relic 
density. Allowing for 3a uncertainty on the WMAP DM abundance, the lower edges of these 
ranges are given to a good approximation by the value obtained without scalar interactions: 
0.51 TeV, 1.7 TeV, 4.4 TeV, 7.6 TeV (0.51 TeV, 2.2 TeV, 9.0 TeV, 21.4 TeV), without (with) 
Sommerfeld corrections for n = 2 complex and n = 3,5,7 real multiplets respectively. The 
upper bound lies from 16 TeV to 60 TeV depending on the model and the perturbativity 
condition one assumes. For a complex multiplet with n = 3, 5, 7 all these values have to be 
reduced by a v2 factor. 

These models are quite predictive. For the inert doublet model, since there are three 
relevant quartic couplings, for a fixed value of the DM mass, there is a two dimensional space 
of values of the quartic couplings that reproduces the observed relic density (see section l3".2.2p . 
Its shape is close to the one of an ellipsoid. This leads to an upper bound on each quartic 
coupling and therefore on the inert Higgs doublet component mass splittings, given in Fig. [2 
As these upper bounds have an asymptotic behaviour for large DM masses, there exists an 
absolute maximum upper bound on each one, from 13 GeV to 17 GeV depending on the 
splitting considered. 

For the multiplets of dimension n > 3, as shown in section 12.21 there is only one relevant 
quartic coupling (the HnH n HlH\ coupling). Therefore the value of this quartic coupling is 
fixed by the mass of the DM, see Fig. [U No mass splitting between the multiplet components 
greater than \{n — l)/2] x 166 MeV are generated in this case. 

For the doublet case the direct detection elastic cross sections can be enhanced by up 

7 for simplicity we assume in the numerical example that N3 is heavier and has little effect on leptogenesis 
8 Note finally that in this framework the N's produce not only a L asymmetry but also a inert doublet 
asymmetry, but for A5 satisfying the direct detection lower bound above this asymmetry rapidly is washed 
away by the A5 driven H 2 H 2 <-► H 1 H 1 processes together with pure Higgs boson self interactions. 
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to 2 orders of magnitude by the scalar coupling contribution, see Fig. The cross section 
is predicted to lie within the range [0.1 — 9.0] • 10~ 9 pb. Consequently they can exceed the 
sensitivity of future planned experiments by a similar amount. For higher multiplets the cross 
section is completely fixed by the DM mass, leading to values which can be ruled out by future 
experiments such as XenonlT, from 0.9- 10 -9 pb to 40 -10 -9 pb. Indirect detection signals can 
also be enhanced by the quartic coupling contributions. For a standard Navarro-Frenk- White 
halo profile, and without any boost factor, the total gamma ray flux is within reach of future 
telescope such as FERMI-LAT. Search for high energy neutrinos from the galactic center is 
complementary to the gamma ray signal. However, without boost, the predicted flux is 1-2 
orders of magnitude below the projected sensitivity of km? size detectors. The antiproton and 
positron fluxes are 3-4 orders of magnitudes below the expected background, see Fig. [JJ Very 
large boost factors would be therefore necessary to have a signal exceeding this background. 
Since the scalar DM models are viable over large ranges of masses, it is clear that some values 
of the mass within these ranges will lie on the top of a Sommerfeld resonance, possibly leading 
to a large boost for all indirect detection signals. A precise determination of these resonances 
potentially relevant for PAMELA and ATIC experiment results is beyond the scope of this 
paper. 

The values of DM masses we have obtained for the higher multiplet are clearly too high 
to allow DM particle production at the LHC collider. It would however be very interesting to 
analyze the possibilities to produce the inert doublet components (in particular the charged 
ones) in the range ~ 0.5 — 1.5 TeV. 

Finally if one adds right-handed neutrinos to the doublet model it is possible to success- 
fully generate in a simple way the neutrino masses, the baryon asymmetry of the universe 
(via leptogenesis in a non-resonant way) and the dark matter relic density, all this in a related 
way and at a scale as low as TeV. Conversely this means that if a second Higgs doublet is 
added to the ordinary type-I seesaw model one can lower the right-handed neutrino masses 
and therefore the leptogenesis scale from 10 9 GeV down to TeV, and at the same time ex- 
plain the observed DM relic density with a TeV scale scalar candidate from the second Higgs 
doublet. 
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A Higher SU(2) L multiplets 

A.l Generators, real and complex multiplets 

In this section, we give explicit expressions for the SU{2) generators in the representation n, 
and define real and complex multiplets. 

Let us consider a multiplet H n in the representation n of SU(2)l. The generators must 
be chosen and normalized so as to satisfy the commutation relations of SU(2), [7o^,Ti ] = 

(n) 

itabcTc i but their form is not uniquely determined. The spherical basis, where the third 
generator is diagonal 

rf* = diag(j n , j n - 1, . . . , -(j n - 1), -j n )) , (68) 

with j n = (n — l)/2, is particularly convenient as the components of the multiplets are 
eigenstates of the electromagnetic charge. The first two generators can be constructed from 
the ladder operators t±^ = ± ir^ which are given by 



' -[(jn-k^jn + k + l)} 1 / 2 ^) , k>0 
[(jn-mn + k + l)] 1 / 2 ^) , k<0 



and = f T +^l > where |ej^) (with k = —j n , —jn + 1, jn) are the basis- vectors. As an 
example, for the triplet case, 

rf > = -*= (-1 lj = | -i -i | t^> = | | . (69) 

These generators are related to the generators (t^ 3 ^) 6c = —ie abc in the cartesian basis by the 
unitary transformation matrix 
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U = ^\i -i | . (70) 



From a group theory point of view, any representation of SU(2) is real, in the sense that 
it is equivalent to its complex conjugate. For the representation n of SU(2), the matrix T n 
which realizes this equivalence, 

T n r^T~ l = -r<»> (71) 

is given in the spherical basis by 

(_!)n+l| e W) j k >0 



T le (n)N 



H) , k<0 



Therefore, for a multiplet H n , the conjugate multiplet H n = T n H* also transforms as the 
representation n under SU(2). 
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We can distinguish between real multiplets for which H n = H n and complex ones for which 
H n ^ H n . A complex multiplet of dimension n contains twice as many degrees of freedom 
as a real multiplet of same dimension. In the spherical basis given above, the components of 
charge Q and — Q of real multiplets are complex conjugate of one another, in the cartesian 
basis, all components are real fields. Finally, it is to check that HnT a H n vanishes for real 
multiplets. 



A. 2 Complex multiplets: mass spectrum and comparison to the real case 



In the section [2~2| we presented in details the mass-spectrum and properties of real multiplets. 
In this section we analyze whether a complex multiplet can still be considered as a minimal 
candidate of dark matter, and how its phenomenology differs from the real case. In the 
spherical basis where the third generator is diagonal, the complex multiplet with Y = can 
be cast in the following form 

/ A C*») \ 



H„ 



A(°) 



(72) 



(-in) / 
2 / 



where the upper index of a component corresponds to its electric charge. Notice that the 
number of independent fields has doubled compared to the real case (which explains the 
absence of a normalization factor l/y/2 in Eq. (|72p ). 

In the case of a complex H n with Y = 0, the lagrangian and potentials given in Eqs. dMS]) 
are not the most general ones anymore, because of the possibility of mixed products between 
H n and H n (like (HnH n ) or (H n T^ H n )). However, the complex multiplet can be decomposed 
into two real multiplets. Indeed, if we define 

1 



A n = —(H n + H n ) 
B n = —j=(H n — H n ) 



(73) 



it is easy to check that A n and B n are real multiplets, that is A n = A n and B n = B n . 
Therefore the most general model with a complex multiplet H n with vanishing hypercharge 
is equivalent to a model with two interacting real multiplets A n and B n . Following the 
minimality criterium that enables us to make a systematic study, we will not pursue the 
details of such a model. 

There is however one case where the decomposition of Eq. (]73p is forbidden, namely when 
H n is charged under some additional gauge group U(1)q'. In that case, the potential of 
Eqs. (JSj) is still the most general one and the A4 and A5 terms do not vanish. In particular, 
A5 generates a mass splitting at treedevel between the components of H n , 



m 



!(A (Q)) 



» 2 + 



^0 



+ {-i) a Q 



(74) 



with Q = 1, . . . , j n and a = 1, 2. As a consequence, half of the charged fields of the multiplet 
are lighter than the neutral component at treedevel, the latter cannot therefore be a DM 



32 



candidate, which rules out the model. At one-loop however, one has to take into account the 
additional splitting generated by the coupling to gauge bosons 

m (aW } ) - m (a(°A = Q 2 AM g , (75) 

with AM g ~ (166 ± 1) MeV. Therefore, the neutral component stays the lightest particle as 
long as 

A5 * 2-2 x 10- (^) , (76) 

where mo = m (A^ ^). In this calculation, scalar sector induced 1-loop corrections have been 
neglected for the following reason. Only A5 induces charge-dependent 1-loop corrections to the 
scalar mass (via loops of the SM Higgs) in the scalar sector so all corrections are proportional 
to A5, A| or A3A5. In the first two cases, it is clearly impossible for 1-loop corrections to 
compensate tree-level splittings because they are proportional to the same coupling. In the 
last case also, A3 would have to be taken far beyond the perturbative regime for 1-loop 
corrections to be non-negligible. 

The constraint Eq. (fT6j) puts a strong upper bound on the only coupling that could impact 
the DM phenomenology of the complex models compared to the real cases since A4 doesn't 
introduce new couplings to the SM particles. As a result, the complex cases are mostly 
equivalent to the real cases, except for the fact that the number of degrees of freedom has 
doubled. As a consequence, the total cross-section of (co) annihilation increases by a factor 
2. As f^DM oc m Q 2 ( see section [3]), for a given relic density, the mass of the DM candidate 
is smaller by a factor \/2 in the case of a complex multiplet compared to the corresponding 
real case. 

B Feynman Diagrams for Annihilation and Coannihilation in 
Inert Multiplet Models 

B.l Inert Doublet Model 

We represent below the annihilation and coannihilation processes. They are organized by 
types of output particles. 
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B.2 Inert Multiplet Model 



Like for the doublet, we represent below the annihilation and coannihilation processes for the 
inert multiplet model of any dimension n. The contributions are organized by type of output 
particles: Figs. [13|14|, and [15] represent the channels involving Higgs particles, the fermionic 
channels and the pure gauge channels respectively. 

One diagram encloses several cases corresponding to all the possible values of Q which 
stands for the absolute value of the charge of the A ±( ^ n-uplet component. Remember that 
Q = 0, 1, . . . ,j n , with j n = (n — l)/2. Notice that some diagrams do not exist for all charges 
or all models. For example, those involving A^ -1 cannot be applied to the Q = case and 
those involving A^ +2 do not exist in the triplet case. Moreover, some interactions are not 
possible because of the absence of coupling. This is the case of e.g. the sixth diagram of 
Fig. [15] for Q = because A does not couple to Z or the photon. 
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Figure 10: Pure gauge (co) annihilation channels 
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Figure 12: (Co)Annihilation channels into fermions for the doublet model. 
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Figure 14: (Co)Annihilation channels into fermions for the Multiplet model. 
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Figure 15: Pure gauge (co) annihilation channels for the multiplet model. 
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